Spinless Fermionic ladders in a magnetic field: Phase Diagram 
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The system of interacting spinless fermions hopping on a two-leg ladder exhibits a series of quan- 
tum phase transitions when subjected to an external magnetic field. At half filling, these are either 
U(l) Gaussian phase transitions between two phases with distinct types of long-range order or 
Berezinskii-Kosterlitz-Thouless transitions between ordered and gapless phases. 



PACS numbers: 71.10.Pm; 71.30.+h 
I. INTRODUCTION 

The behavior of interacting electron systems under the 
action of an external magnetic field, affecting orbital mo- 
tion of the particles, is a subject of intense research of the 
last several decades. It has long been established that 
the magnetic field has a dramatic effect on all properties 
of the system. Even in the absence of interaction, the 
spectrum of the free Fermi gas is modified and exhibits 
Landau quantization^ in the continuum or the Hofstadtcr 
spectrum^ on a two-dimensional (2D) lattice. The most 
spectacular consequence of this phenomenon is the Quan- 
tum Hall Effect (QHE). 3 Electron-electron interaction 
compounds the complexity of the problem, giving rise 
to the Fractional Quantum Hall Effect (FQHE)ii 

The current theoretical understanding of the effect of 
the magnetic field on the properties of electron systems 
was achieved by a combination of various methods and 
techniques, each of which is strictly speaking only appli- 
cable in a certain parameter range. What is still lacking 
is a comrehensive approach that would unify all of the 
different aspects of the problem into a single coherent 
picture. Perhaps at present such a goal is too ambitious. 
However, as a first small step in this direction, one can 
ask whether such a comprehensive approach can be for- 
mulated for a simpler model, that would on one hand be 
a problem of interacting electrons in the magnetic field 
and on the other hand would, at least in principle, allow 
a generalization towards the original problem. 

The simplest model of interacting fermions that in- 
corporates orbital effects of an external magnetic field 
is that of spinless fermions hopping on a two-leg ladder. 
This model is simple enough not to exhibit the multitude 
of small gaps in the single-particle spectrum (character- 
istic of the Hofstadter problem 2 ). Yet the magnetic flux 
piercing the plaquettes of the ladder changes the ground 
state properties of the system giving rise to non-trivial 
phases and inducing quantum phase transitions. 

Ladder models^* 6 - occupy a special place in the field of 
strongly correlated electron systems. On one hand, they 
describe (at least within some range of temperatures) 



behavior of many naturally found compounds 7 , includ- 
ing carbon nanotubespiS as well as artificially manufac- 
tured structures^ On the other hand, they provide a 
fertile ground for application of theoretical techniques 
developed for one-dimensional systems^ - i.e. non- 
perturbative approaches leading to asymptotically exact 
results. Furthermore, they are a first step for various 
attempts at generalization of the lore of physics in one 
spatial dimension to higher-dimensional problems. 12 

As is often the case, situations where the number of 
particles is commensurate with the lattice attract the 
most attention since only then can long range order 
(LRO) develop. We have previously showni^ that at a 
quarter filling, the presence of the flux results in exciting 
effects which do not exist in the absence of the flux. In 
particular the uniform external magnetic field can lead 
to a staggered flux (or orbital antiferromagnet) phase, 
which furthermore has fractionally charged excitations. 

The purpose of the present paper is to describe the 
full phase diagram of interacting spinless fermions on the 
two- leg ladder at 1/2-filling in the presence of an ex- 
ternal magnetic field. To drive a system to criticality 
by applying the magnetic field is an intriguing possibil- 
ity which should be more accessible in experiment than 
varying coupling constants. The phase diagram displays 
a multitude of quantum phase transitions induced by the 
flux. There are two types of these phase transitions: (i) 
Berezinskii-Kosterlitz-Thouless (BKT) transitions^ be- 
tween ordered and disordered ground states, and (ii) 
U(l) Gaussian phase transitions between different or- 
dered ground states. Here we choose the fermions to 
be spinless in order to eliminate Zeeman splitting and 
focus on the orbital effects of the magnetic field. The 
case of spin- 1/2 particles will be discussed in a separate 
publication. 

Traditionally 15 - 16 ladder models have been treated in 
two complimentary approaches. On one hand, one can 
start with the model of two decoupled chains, define the 
low energy effective theory for each of them, and then 
treat both the single-particle transverse hopping and two- 
particle inter-chain correlations perturbativelyJ- 7 . On the 
other hand, one could start with the exact single particle 
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basis (given by two bands) and then proceed with the 
corresponding low energy limit & For the major part of 
this paper we will be using the latter approach which al- 
lows us to treat intra- and inter-chain processes on equal 
footing. However, it is well known- (at least in the ab- 
sence of the magnetic field) that for weak enough inter- 
chain tunneling there exists the phenomenon of Ander- 
son confinement, i.e. the suppression of the inter-chain 
single-particle tunneling by intra-chain two-particle cor- 
relations. This effect can be seen in either picture, but it 
is more intuitive to discuss it within the chain approach. 
Both ways should produce the same physical results, but 
understanding the relation between the two approaches 
allows us to establish the limits of applicability of the 
effective low energy theories that can be derived within 
either picture. Moreover, we discuss how the Anderson 
confinement regime is affected by the magnetic flux. 

The remainder of the paper is organized as follows. 
We start by defining the microscopic Hamiltonian of the 
model and proceed directly to the results, discussing the 
phase diagram and the other physical properties of the 
model. Then we outline the details of the calculations 
within the weak coupling (bosonization) approach. In 
Section IIHI we derive the effective low energy theory in 
the band picture. In Section IIVI we turn to the chain 
picture and discuss its relation to the band approach. 
Section is devoted to the strong coupling limit of our 
model and is followed by a brief summary of the results. 
Mathematical details are relegated to Appendices. 



models of spinless fermions and Hubbard-like chains of 
spin-1/2 particles. In our case however, the SU(2) sym- 
metry is explicitly broken in the Hamiltonian Q by the 
inter-chain hopping and the V interaction term (the for- 
mer is analogous to a Zeeman energy due to an external 
magnetic field alng the x-axis, while the latter is a coun- 
terpart of an exchange anisotropy along the z-axis). 

The external magnetic field B is introduced by means 
of the Peierls substitution^ In the Landau gauged with 
the vector potential A = B(—y, 0, 0) the transverse hop- 
ping term is independent of the field, while the longitu- 
dinal hopping amplitude can be written as 



(2) 



where / is the magnetic flux through the elementary pla- 
quette in units of flux quantum </>o = hc/e. Expressed in 
terms of the flux the model is explicitly gauge invariant. 
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II. MODEL AND RESULTS 

In this section we present our results. We start by 
defining the microscopic Hamiltonian of our model and 
proceed to discuss the zero-temperature phase diagram. 



FIG. 1: Cartoon depictions of possible ordered phases at 1/2 
filling. The dots represent excess fermion dimerization on the 
sites. The ellipses represent excess occupation on the bonds. 
The arrows represent local currents. Note that the OAF and 
the BDW coexist - see Section UTTl for an explanation. 



A. The Hamiltonian 

We consider a tight-binding model of spinless fermions 
on a two-leg ladder described by the Hamiltonian 

H= - ~^2[t\\(a)cl a c i+ha + h.c. - t x ^2 c la°i -o 



(1) 



Here Ci t „ is the electron annihilation operator on the 
chain a = ± at the site i; nj CT = c] a c%. a are the occupa- 
tion number operators; t± and tn are the transverse and 
longitudinal hopping amplitudes, respectively. The last 
two terms in Eq. describe nearest neighbour inter- and 
intra-chain interactions. The chosen form of short-range 
interaction is quite representative because it reflects the 
generic symmetry of the ladder and yields the most gen- 
eral effective field theory in the low-energy limit. Our 
notation reflects the well-known analogy between ladder 



B. Phase Diagram 

At half filling (one fermion per rung) the model Q is 
characterized by a rather rich phase diagram. Depending 
on the values of microscopic parameters the ground state 
of the model may posess true long range order. Possi- 
ble ordered phases are illustrated in Fig. The cartoons 
show a strong-coupling picture of the phases: charge den- 
sity waves (CDW), where particles are localized on sites 
of the ladder; bond density waves (BDW) with dimerized 
links along the chains; and the orbital antiferromagnet 17 
(OAF), sometimes referred to as the staggered flux phase 
or a d-density wave, where the particle density remains 
uniform, but there exist non- vanishing local currents that 
have opposite directions on alternate bonds. Notice that, 
in the spin language, the Relative CDW is similar to a 
spin density wave (SDW) with spins polarized along z (a 
Neel state SDW Z ), whereas the OAF is equivalent to a 
SDW y . 
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In addition, the model (JTJ allows for various phases 
that do not posess long-range order. Using the afore- 
mentioned analogy with the Hubbard chain, we may dis- 
cuss the model Q in terms of "spin" (or "relative" ) and 
"charge" (or "total") sectors. In Section Till Bl we show 
that, in the low-energy limit, the "charge" and "relative" 
sectors of the model asymptotically decouple. In all of 
the ordered phases both sectors are gapped. However, it 
is possible to have a phase where only one of the sectors 
acquires a spectral gap. Phases where only the "charge" 
sector is gapped, irrespective of the type of dominant 
correlations, we will call the Mott Insulator (MI). The 
cases where the gap exists in the "relative" sector only 
will be called the Luther-Emery Liquid (LEL). 21 Finally, 
when both sectors are gapless, the system represents a 
Luttinger Liquid (LL). 
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FIG. 2: Phase diagram at B = (after Ref. EH- Phase 
boundaries correspond to the lines V = 0, U/2V = —2 + r 2 
and U/2V = ~r 2 . 

The phase diagram in the absense of the flux is 
knownii£i±& For the sake of clarity we include it in Fig. 
This phase diagram is valid for sufficntly large values 
of t± where derealization of the fermions across the 
rungs suppresses the CDW phase (which happens in 
the absence of the inter-chain tunneling). There are 
two ordered phases: (i) for purely repulsive interactions 
U, V > 0, one has a relative CDW as to be expected (see 
Section 0) ; (ii) for repulsive interchain interaction and 
not too strong attractive in-chain interaction, the ground 
state is the orbital antiferromagnet, 1 ^ 22 The phase di- 
agram in Fig. |21 was obtained within a weak-coupling 
bosonization approach. The phases do exist when the 
coupling becomes strong, however the exact location of 
the phase boundaries might change. 

Once the magnetic field is applied, the system may 
exhibit additional phase transitions. In Fig. [3] we plot 
the entire weak-coupling phase diagram for the model 
at half-filling and sufficiently large t± (see next subsec- 
tion). The magnetic flux varies along the vertical axis, 
so that the diagram in Fig. [2] corresponds to the bottom 
axis of Fig. [3J The ratio of the microscopic interaction 
parameters of the Hamiltonian Q, U/2V, is represented 
in Fig. 13 through the angular variable 6 — tan -1 U/2V 



which varies along the horizontal axis. Four different re- 
gions corresponding to various signs of the constants U 
and V are indicated on the upper horizontal line in Fig.|3 
The analytic description of the phase boundaries, based 
on the weak coupling theory, is given in Appendix ^ 
The position of the boundaries depends on the applied 
field, the ratio U /2V (i.e. ff) and the ratio of the hopping 
parameters, r = t±/to- We plot the phase diagram as a 
function of 9 and flux for r = 0.25, the value is arbitrary 
but representative as long as r exceeds a possible gap 
in the "relative" sector. At other values of r the topol- 
ogy of the phase diagram and classification of the phases 
do not change qualitatively. Similarly, if we modify our 
model ^ to include other short-range interaction terms, 
the only effect on the phase diagram would again be just 
the shift of the phase boundaries. 

Let us now describe the phase transitions induced by 
the applied field (i.e. the vertical direction in Fig. [3}. 
We assume that the inter-chain hopping parameter r is 
not too small (see next subsection). Note, that since 
the model Q is invariant under the transformation / — > 
1 — /, a — > — cr,we only need to consider the flux within 
the range < / < 1/2. Moreover, when the flux is large 
enough, sin 2 nf > I — t 2 , there is a band gap in the single 
particle spectrum of the model. That state is largely un- 
affected by interaction effects (at least within the weak- 
coupling limit), and thus we will restrict our discussion 
to smaller values of /. While the weak-coupling approach 
can not be trusted at fields too close to the band gap limit 
(since the Fermi velocity becomes too small) , we continue 
the phase boundaries up to that point. All phase transi- 
tions of interest happen sufficiently far from that region. 
Apart from a brief discussion in Sec. we will not con- 
sider the details of the transition to the band insulator 
in this paper. 

The most interesting features of the phase diagram in 
Fig. are a sequence of U(l) phase transitions between 
different ordered states and reentrant transitions. Un- 
derstanding of these transitions is based on the fact that, 
as shown in Section 1111 HI in the low-energy limit the 
"charge" and "relative" degrees of freedom of the model 
decouple, and each sector is described by a sine-Gordon 
model (see Eq. (|I2|) and Appendix 1X1). Phases with LRO 
correspond to strong-coupling regimes in both sectors. 
The phases whose order parameters are mapped onto 
each other under a sign change of the corresponding cou- 
pling constant (the amplitude of the cosine term) are mu- 
tually dual. The associated U(l) Gaussian criticality oc- 
curs at the self-dual lines, i.e. when the one of those cou- 
pling constants vanish. Such a duality is commonplace 
in low-energy effective theories, indeed more complicated 
non-Abelian dualities were found recently in the SU(4) 
Hubbard modeL 2 - However, it has also been recently 
shown 24 that for certain ladder models the (Abelian) du- 
ality between different phases turns out to be not only a 
symmetry emerging in the low-energy limit but an exact 
property of the underlying microscopic model. 

Gaussian transitions occur in two domains of the phase 
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FIG. 3: The weak-coupling phase diagram in the magnetic field. We plot the flux along the vertical axis and the angle 8 
(defined as 6 = tan^ 1 U/2V) along the horizontal axis. As the phases depend on the signs of the interaction parameters, 
they are indicated at the top of the diagram. Ordered phases are illustrated pictorially in Fig. Q The corresponding order 
parameters are listed in Table The disordered phases are characterized by dominant correlations (indicated in parentheses). 
For large values of the flux (sin 2 irf > 1 — r 2 ), there is a band gap in the non-interacting picture. The thick solid lines (blue and 
green online) represent U(l) Gaussian transitions between mutually dual ground states with long-range order, and the thick 
dotted lines (black and red online) are Berezinski-Kosterlitz-Thouless phase transitions corresponding to opening of a gap in 
one of the sectors. 



diagram: (i) for repulsive intrachain interaction (V > 0) 
and weak attractive interchain interaction (U < 0), and 
(ii) for U > and weak V < 0. In the first case (i), at zero 
flux the system is a MI. The dominant (longest-range) 
correlation turns out to be that of the 2k p component 
of the total charge density (hence the label "(CDW)" in 
Fig. |21 see Section 1111 PI for details on dominant corre- 
lations). As the flux is increased, the system becomes a 
relative BDW (via the BKT transition where a gap opens 
in the "relative" sector of the effective theory). Further 
increase of / drives the system through the t/(l) transi- 
tion to a CDW phase. At larger values of the flux (ap- 
proaching the band gap limit) values of flux, the system 
eventually becomes an OAF, again through the Gaussian 
transition. 

The second transition (ii) occurs at small values of the 
flux. The zero-field ground state is an OAF. As we turn 
on the flux, the system undergoes a U(X) transition to- 
wards a relative CDW state. Further increase of the flux 
results in closing of the gap in the relative sector (via the 



BKT transition) and the system becomes a MI, but now 
the dominant correlation is that of the 2kp component 
of the transverse bond density, labeled by "(BDW)". 

This latter transition turns out to be reentrant. As the 
flux is further increased the system returns (again via the 
BKT transition) back to the relative CDW state. There 
is another example of a reentrant transition in the phase 
diagram, for U < and small V the zero-field ground 
state is a Luttinger Liquid, which once subjected to the 
external field first becomes a LEL by opening a gap in 
the "relative" sector, and then at higher field comes back 
to the LL state in which the most singular fluctuations 
are those of the pairing operator at momentum tt — 2kp. 
In the LEL phase the "relative" sector is gapped and the 
"charge" sector is characterized by the dominant corre- 
lation of the pairing operator at zero momentum. 

Finally there is a large part of the phase diagram which 
is robust against the application of the external field. 
When both inter- and intra-chain interactions are at- 
tractive, the LL (that is the zero-field ground state) is 
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mostly unaffected by the field. More interesting is the 
situation when both interactions are repulsive. The zero- 
field ground state is the relative CDW. It turns out that 
this long-range order survives under the application of 
the field (except possibly for the transition to the MI for 
weak V discussed above). 



C. Commensurate-Incommensurate transition 

The above phase diagram breaks down if the param- 
eter h — [sin 2 (7r/) + r 2 ] 1 / 2 , which in the noninteract- 
ing case determines the splitting between the Fermi mo- 
menta of different bands, is too small (see Sections IIIII 
and IIV(I . Then, the part of the phase diagram that cor- 
responds to attractive inter-chain interaction U < ac- 
quires additional ordered phases. This is the result of 
additional inter-band scattering processes that at larger 
h violate momentum conservation in the low energy ef- 
fective theory (based on the two-band description). The 
latter issue reveals the dichotomy between the two start- 
ing points already mentioned in the Introduction: chain 
basis versus band basis. If one starts with a solution 
(however complete) for two independent chains and then 
tries to take into account the inter-chain hopping (as 
well as the flux) in perturbation theory, then the pro- 
cesses mentioned above are present in the theory from 
the beginning. In the case when these processes gen- 
erate a gap in the spectrum of relative degrees of free- 
dom, a finite splitting of the Fermi momenta would not 
take place unless the parameter h exceeds its critical 
value comparable with the gap. This is the well-known 
commensurate-incommensurate transition^ As the pa- 
rameter h increases further, the "two-chain" approach 
fails because renormalization of the parameters of the 
theory becomes sizable at sufficiently large h. At that 
point one would be forced to start with the exact, two- 
band single-particle spectrum of the ladder. However, 
this would seemingly neglect the processes in question as 
they appear to violate momentum conservation. 

In Section Hvl we discuss the relation between the two 
approaches to ladder problems and show that if one uses 
either approach properly, then the final result is inde- 
pendent of the starting point, as should be expected. 
The new phases at U < naturally emerge through the 
commensurate-incommensurate transition^ 

We shall also show that, regardless of the starting 
point, some properties of the system are not accessible 
within the effective low energy theory. The quantity in 
question is the diamagnetic (or persistent) current, which 
turns out not to be an infra-red property. All electrons 
participate in this current. In particular, the curvature of 
the single-particle spectrum at the Fermi points becomes 
important, so that linearization of the spectrum, being 
the usual prerequisite in the derivation of any effective 
low-energy theory, completely destroys this effect. Con- 
sequently, within the bosonization approach in the con- 
text of the fermion ladder, it is impossible to describe the 



analog of the Meissner effect that can be seen in bosonic 
ladders^ Details are presented in section HVI 



III. LOW ENERGY EFFECTIVE THEORY: 
BAND BASIS 

In this Section we derive the effective low energy theory 
for the model taking the exact single particle spec- 
trum as our starting point. As mentioned above, there 
exists an alternative approach, which starts with discon- 
nected (but interacting) chains. The relation between 
the two will be discussed in the next Section. 



I II 




FIG. 4: Possible types of the single-particle spectrum as a 
function of in-chain momentum. For any given r, increasing 
the flux will eventually open a band-gap in the non-interacting 
spectrum. Close to this transition one of the two bands is 
almost empty while the other is almost full. At this point, 
curvature effects of the spectrum become important; these 
are beyond the scope of this paper. 



A. Single-particle spectrum 

The single-particle part of the Hamiltonian can be 
diagonalized by the unitary transformation 

ci(fc) = u k a k + v k (3 k , 

ci{k) = v k a k - u k f3 k (3) 

where the "coherence factors" u k , v k [which are positive 
as the signs are written explicitly in Eq. (J2J)] are given 
by (if not stated otherwise, in this section we will mea- 
sure the momentum k in units of the inverse longitudinal 
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lattice spacing, l/a) 



absence of interaction) is given by 



sin k sin irf 



V sin 2 k sin 2 irf + r 2 

sin k sin irf 
y sin 2 A; sin 2 7r/ + r 2 



(4) 



The resulting single-particle Hamiltonian describes two 
bands 



Ha = ^ [ e a a I a fe + e pPkh 

k 



with the spectrum 



£q(/3) 



= -2tn 



cos k cos 7r/ ± y sin fc sin nf + t 2 



(5) 



(j're 



dk 

-2t sinnf / — < cos fc [n a (&) + n,g(fc)] 

Z7T 



! k cos nf [n a {k) — ng(fc)] 
\J sin 2 A; sin 2 7r/ + t 2 



(8) 



where n a ^(k) — (c^ /3 - ) (fc)c Q ,( / g)(A;)) are the occupation 
numbers of the two bands. The current © is a periodic 
function of the flux with a period Af = 1. In Fig [S] we 
plot the flux dependence of (j re i) within a single period. 
Notice that the current changes its sign under transfor- 
mation / — > n — /; at / = 1/2 (jrci) = due to the 
recovery of time reversal symmetry at this point. 



In the absense of the flux the coherence factors are inde- 
pendent of momentum (it 2 = v 2 = 1/2) and the Hamil- 
tonian Hq consists of the usual symmetric and antisym- 
metric bands each with the cosine spectrum, split by 2t± . 

In the presence of the flux the spectrum can take one of 
the four typical shapes depending on the value of r and 
the flux. These are illustrated in Fig.0] For completeness 
we include a detailed discussion of the properties of the 
spectrum as a function of flux and interchain hopping in 
Appendix [5] 

The half-filled ladder is characterized by zero chem- 
ical potential. When the spectrum © takes the form 
depicted in the two bottom graphs in Fig. 0] the system 
at fi = is a band insulator. In that case interaction 
effects (as well as the external field) are not expected to 
drastically change the nature of the ground state of the 
non-interacting system. We will not discuss that case in 
the present paper. 

The top two graphs in Fig. 0] describe the "metallic" 
phase of the non-interacting system. In this case at half 
filling both bands are partially filled and each band is 
characterized by its own Fermi momentum k F ^ satisfy- 
ing kp + kp = n. In what follows we will use the notation 
kp = kp (so that k F =n — kp) with 



cos kp = y sin 2 nf + t 2 . (6) 

In the presence of the magnetic flux there exists a fi- 
nite diamagnetic current in the ground state of the sys- 
tem. The current operator along the oriented link be- 
tween sites n and n + 1 of the chain a is 

j n ,a = -ito (e l7rfa c} h<T c n+1 . a - h.c.) . (7) 

This current flows in opposite directions on the two legs 
of the ladder, so that the total current jtot — jn,+ + jn,- 
will have zero expectation value, while the expectation 
value of the relative current jrei = jn,+ — in.- (in the 
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FIG. 5: Diamagnetic current as a function of flux in the ab- 
sence of interaction. Only one period in / is shown. The 
cusps correspond to the band gap opening. 

In the limit r — > 0, when the ladder decouples into two 
completely disconnected chains, the appearance of the 
flux in the Hamiltonian is a gauge artifact. Indeed, 
a careful evaluation of the integral in Eq. @ will show 
that (jrei) — at t = 0. Expanding Eq. (JHJ for small 
/,r <C f and recovering the dependence on the lattice 
spacing a, one finds that 

(> ei ) = 2/r 2 [f + 0(/ 2 ,r 2 )]. (9) 

Despite being small in this limit, the diamagnetic cur- 
rent is not an infra-red phenomenon. Its dependence on r 
is apparently the effect of a finite curvature of the single- 
particle spectrum. Notice, that as seen from Fig El Orel) 
is non-zero even in the insulating phase. Thus the current 
is a non-universal quantity contributed by all electrons 
and not only those in the vicinity of the Fermi points. 
Consequently, effects related to such a persistent current 
can not be addressed in terms of any Lorentz-invariant 
effective low-energy theory (we will further comment on 
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this issue in Section [lV|l . Thus, at present we are un- 
able to calculate the effect of the interaction on the dia- 
magnetic current. However, it is clear that even in the 
presence of interaction the current will still persist and 
all the correlation-related phenomena discussed in this 
paper will coexist with it. 

B. Interaction Hamiltonian 

Now we are going to apply the standard rules of 
Abelian bosonizationi£*i& to derive the effective low- 
energy theory. First, we will assume that the Fermi en- 
ergy is sufficiently far from the bottom of the /3-band. 
Then we linearize the two-band spectrum Eq. (JSJ in 
the vicinity of the four Fermi points, ±k F and ±A;£. 
The associated low-energy degrees of freedom are de- 



scribed in terms of smoothly varying chiral (right and 
left) fermionic fields, R a r^(x) and L a ^(x). This de- 
fines the continuum limit of the model in which the 
non-interacting part of the lattice Hamiltonian, including 
both the inter-chain hopping and the coupling to the flux, 
transforms to the kinetic energy of the chiral particles: 




where vp — 2t^a sin k F / cos nf is the Fermi velocity 
which at half filling is the same for both bands. 

Specializing to the vicinity of the four Fermi points 
in the coherence factors Eqs. (@J we find the low-energy 
correspondence between the original lattice operators Ci jCr 
and the chiral fields R v and L v . Then, the interaction 
terms in the model become 



Hint ~ a ^ |gi (: JruJlo ■ + ■ JrpJlp '■) + 92 ('■ JruJr/3 ■ + ■ JloJl/3 '■) + 93 ('■ JruJl/3 ■ + : JloJrp '■) 
+.g 4 (: RiL a R\L (i : + : L^R^R^ :) - g 5 (: R) a l) a RpLp : + : R a L a R^L^ 



+96 {■ R^ a {xi)R^ a {x l+1 )Lp{x i )L l3 {x l+ i) : + : 

+36 L i ' a {xi)L^ a {xi + i)Rf 3 {xi)R (j {x l+ i) : + : L a {xi)L a {xi + i)^(xi)R^{x i+ i) j. (10) 

where Jr„ =: R\,R V : and Jl„ =: L\,L V : are the chiral densities of the right- and left-moving fermions with the band 
index v (the symbol stands for normal ordering). 

The first three terms in Eq. QlOfl. characterized by coupling constants 31,32 and 33, describe the density-density 
interaction, whereas terms with amplitudes 34, 35 and ge correspond to the interchain umklapp, interchain back- 
scattering and in-chain umklapp terms, respectively. Explicit expressions for the gi in terms of the original microscopic 
theory is given in Appendix 1X1 The coupling constants depend on the interaction constants of the microscopic model 
JTJ and, through the coherence factors Eq. (J3J, on the external flux. The latter dependence plays an important role 
because it is responsible for the sequence of phase transitions, described in Section Til Bf) . that are not accessible at 

/ = o. 

The interaction Hamiltonian Hi nt in Eq. (|10|> is the most general form of four-fermion interaction in the band 
representation, consistent with momentum conservation (modulo the reciprocal lattice vector). All other terms contain 
strongly oscillating exponentials and thus do not contribute to the low-energy theory. In particular, this argument 
applies to the term 

g 7 e 2i ^- k >RiL a LlR + h.c. (11) 

Note, that if one starts building the low energy theory approximating the ladder by two uncoupled chains, then the 
Fermi momenta of the two bands are equal and the above term should be included in Eq. i|l(J|) . We will discuss this 
term and the relation between the two approaches to bosonization in ladder models in Section IIVI However it is 



immediately clear, that omitting Eq. (|llf> from Eq. I|1U|) can only be valid at long distances \x\ ^ \k F — k F \ 1 or, 
equivalently, at low energies \uj\ ^ vp\k F — k F \. 

Now we bosonize the theory in the standard manner (our conventions are outlined in Appendix . As usual, the 
density-density terms (represented in Eq. i|10|) by 31, 32, and 33) renormalize the Fermi velocities and the scaling 
dimensions of the vertex operators. Introducing symmetric and antisymmetric combinations of the bosonic fields, 
4>± = (0q i 4>p)/^/2, we diagonalize the quadratic part of the effective bosonized Hamiltonian. The latter is then 
represented by two sine-Gordon models defined in the symmetric and antisymmetric sectors which are coupled by the 
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in-chain umklapp term g§: 



K + {d x O + f + ^-{d x cj> + f 



.94 



2tt 2 c 



2 



■ cos V8ir(f> + 



.95 



2tT 2 < 



8tt0_ 



.9(5 



cos v 87T(/> + cos y/8n9- . (12 



The cosine terms in Eq. ()12|) , when relevant (these are 
the cases K + < 1, X_ > 1, respectively), are responsible 
for a dynamical generation of a mass gap in the corre- 
sponding sector and, therefore, for the U(l) phase tran- 
sitions described in Section III Bl For weak interaction, 
ISil/Vup <C 1, the "Luttinger liquid" parameters K± are 
close to unity (see Appendix ^J. Consequently, the co- 
sine terms having scaling dimensions 2K + and 2/ K- in 
the symmetric and antisymmetric sectors, respectively, 
are nearly marginal. The g§ term that couples the two 
sectors is therefore strongly irrelevant because of its scal- 
ing dimension 2K+ + 2/K- ~ 4. The only situation when 
the g$ term may become important is the case when one 
of the sectors is gapped, while the amplitude of the co- 
sine term in the other sector vanishes, i.e. either g± = 
or <75 = Q. In this case the ge term can generate the 
missing cosine in the Gaussian sector and even make the 
latter massive. This mechanism was recently discussed 
in Ref . in the context of the Mott instability of a half- 
filled fermionic ladder with U = 0. Since in our model 
the presence of the 174 and g$ terms is generic, and the 
lines 34 = and g§ — characterize the phase bound- 
aries, the only effect of the in-chain (t/g) Umklapp scat- 
tering would be to modify the equations that determine 
the phase boundaries without changing the topology of 
the phase diagram. Being interested in the description 
of distinct phases rather than their precise location, we 
will ignore the g§ term in the remainder of this paper. 

Thus the effective low energy theory for our model 
Eq. m consists of two asymptotically decoupled sectors, 
each being a sine-Gordon model. In the case when a 
strong-coupling regime developes in either sector, a mass 
gap gets generated in the spectrum, and semiclassical so- 
lutions of the equations of motion describe locking of the 
bosonic field in one of the infinitely degenerate minima 
of the cosine potential. Physical quantities evaluated on 
such solutions may either vanish or acquire a nonzero ex- 
pectation value. The former would mean that the quan- 
tity in question is characterized by exponentially decay- 
ing correlations. In contrast, the latter corresponds to 
long-range correlations. Since local operators of the the- 
ory have a multiplicative structure, they can indeed serve 
as order parameters if gaps are generated in both sec- 
tors simultaneously. The multiplicity of the actual values 
that the order parameter would take on the semiclassical 
solutions, differing by a period of the cosine potential, 
determines the degeneracy of the ordered ground state. 
The latter always appears to be associated with a unit 



cell doubling and is The complexity of the formulae, 
relating the four coupling constants in (|12(l to the two 
interaction parameters in the original Hamiltonian 
as well as the magnetic field, leads to a rich phase dia- 
gram, as we will now demonstrate using the just outlined 
strategy. 



C. Ordered Phases 

As there are four distinct Fermi points in our model, 
any local operator will contain four dominant Fourier 
components 

0(x n ) = O (x) + (-l) n O n (x) 

+ cos{2k F x n ) [0 2kF (x) + (-l)"a-2fc F (i)] • 

Here Oq is the smooth part of the operator 0(x n ) cor- 
responding to characteristic momentum q ~ 0; is the 
staggered part contributed by momenta q ~ w, which can 
originate from some inter-band pairing; the components 
02k F and OTr-2k F can be present due to in-band pair- 
ing. At half filling, it is only the staggered part that can 
acquire an expectation value and serve as an order pa- 
rameter; however in some of the gapless phases dominant 
correlations may occur at 2kp or 7r — 2kp rather than ir. 

Local operators of interest in the case of half-filled lad- 
der are listed in Table [I] which includes the microscopic 
lattice definitions and the bosonized form of the domi- 
nant Fourier components. These are given up to mul- 
tiplicative factors; we preserve, however, prefactors pro- 
portional to the magnetic flux to make clear which quan- 
tities do not exist in zero field limit. We also indicate the 
LRO that appears when order parameters (first five op- 
erators in the Table HJ- which also correspond to the five 
"cartoons" in Fig^l acquire nonzero expectation values. 

An interesting observation that can be made from Ta- 
ble|I|is that the OAF and the BDW are both proportional 
to the same low-energy operator, implying that the two 
phases coexist. However, the OAF can exist already at 
zero flux whereas the BDW order parameter is propor- 
tional to the flux (at small /). This coexistence can be 
understood by noticing that at / ^ the BDW order 
parameter is defined in a gauge invariant way, i.e. with 
flux-dependent phase factors explicitly included into its 
definition [we remind that we have chosen the longitudi- 
nal Landau gauge, see Eq. J2J)]. As a result, the BDW 
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TABLE I: Local operaors in the half-filled ladder 



Local Lattice 
Operator Definition 

J± —itx \c\c-2 — ft.cj 

p- c[ci — c\c2 

p\\,+ e l7ri c\{x n )d{x n+1 ) + e~ lw] c\(x n )c 2 {x H+1 ) + h.c. 

p\\- e l7T f c\(x n )c 1 (x n+1 ) - eT mS cl(x n )c 2 (x„ +1 ) + h.c. 

p+ eld + c\c 2 

p± c\c2 + cla 

O sc dC2 



Dominant Bosonized Ordered 

Component Form Phase 

7r cos \/2tt(/> + cos OAF 

7T cos V2^4>+ sin V2kQ- Rel. CDW 

7r sin 7r/ cos \/2TT(f>+ cos \/2tt8- BDW 

tt sin y/2H<j) + sin V^tt0_ Rel. BDW 

7r tan nf sin \/2~TT(f)+ cos \/2n6- CDW 

2kF sin \/2~K(f> + cos \/2~K(f>- 

2k,F cos %/27r(/)+ sin v / 27T0- 
7T - 2fc F ie ,;v ^ fl + cosfv 7 ^- - (tt - 2k F )x] 

tan 7r/e ,V5¥e + cos v^fl- 



operator describing dimerization of the two chains with 
zero relative phase acquires an admixture of the stag- 
gered relative current, proportional to the flux at / <C 1. 
This admixture actually represents the longitudinal part 
of the OAF order parameter (which is identical to J± 
by current conservation). The very appearance of such 
an admixture is a consequence of the explicit breakdown 
of time reversal symmetry, caused by the external flux, 
which is superimposed on the spontaneous breakdown of 
this symmetry in the OAF phase. 

Another ordering in which the flux plays a crucial role 
is the CDW phase. As already mentioned, the ground 
state of the ladder at / = and not too small r does not 
display this type of LRO as the interchain hopping tends 
to prevent double occupancy of the rungs. It is a curious 
fact that under application of the flux this state can be 
recovered due to a similar, although more subtle, admix- 
ture with a staggered flux phase. Indeed, the bosonized 
low-energy projection of the operator p + has the form 
/^diag, where the operator 

Jdiag ~ «(-!)" cr (4 iCT c n+ i _ CT - h.c.) , (13) 

(7 

represents an order parameter for an OAF state with lo- 
cal currents effectively flowing across the diagonals of the 
plaquettesi2£ The scalar nature of the CDW under time 
reversal is not violated for the reason already mentioned 
in the preceeding paragraph. 

Let us now turn to the derivation of the phase diagram. 
For the cosine terms to become relevant and generate a 
gap in the spectrum, the "Luttinger liquid" parameters 
K+ and K- should be smaller and larger than 1, respec- 
tively. According to the definition of K±, Eas. (|A4|l . this 
translates into the following conditions on the parameters 
of the theory: 

5 i + <73>0; 5i- 53 <0. (14) 

Strictly speaking, these conditions are valid only to first 
order in the Kosterlitz-Thouless RG equations. When 



the couplings <?j (i = 1, 3, 4, 5) are all of the same order, 
there are important renormalizations of the parameters 
K± emerging in the second-order— 6 . This means that 
the exact positions of the phase boundaries depend also 
on (?4 and c/5. These corrections, however, do not cause 
qualitative changes in the overall phase diagram and can 
therefore be neglected in the leading order. For this rea- 
son, when drawing conclusions on relevance or irrelevance 
of various perturbations, we will resort to an estimation 
of their Gaussian scaling dimension. 

In the effective Hamiltonian i|12|) there are two cosine 
terms with amplitudes g^ and 55. Both terms have the 
same period in their respective variables that defines the 
values of the fields </>+ and 0_ for any semiclassical solu- 
tion. Depending on the sign of 54 the field 0+ may take 
one of the two possible sets of values, po — or 
(p„ = \j 'it + 'tt /2 (n — 0, ±1, . . .). Similarly, 0- may 
take one of the above values depending on the sign of 35. 
Consequently there are four possible ordered phases. 

(i) If both (74 and 55 are negative, then the semiclassical 
solutions are </>+ = tp„ and 9- = <po- Of all the operators 
listed in Tablc^only the staggered component of the total 

( s) 

charge density p + has a non-zero expectation value on 
the above solution. Therefore, the conditions 54(5) < 
define the charge density wave (CDW). This phase exists 
only in the presence of the magnetic field, which mixes it 
up with an OAF phase, as explained previously. 

(ii) For g4 > and g§ < we find the orbital anti- 
ferromagnet (OAF), since now <fi + = 0- = ipo and the 
staggered component of the inter-chain current J± gains 
the expectation value. In contrast to the quarter-filled 
casein, the OAF phase exists even in the absence of the 
magnetic fieldi^SiSS We will clarify this issue in Sec. [V] 
At / ^ the OAF coexists with BDW, as we already 
mentioned. 

(iii) When both interaction constants are positive the 

(s) 

staggered component of the relative charge density p_ 
has an expectation value (since in this case 4>+ = <po and 
6>_ = (p w ). We call the corresponding ordered phase a rel- 
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ative charge density wave (Relative CDW) . The magnetic 
field has little effect on this phase except for the exact 
location of the phase boundary on the phase diagram, 
which is beyond the scope of this paper. 

(iv) Finally, if 54 < and 95 > then 4> + = 9- = ip^ 
and the staggered component of the relative longitudinal 
bond density p}, _ acquires a non-zero expectation value 
yielding the relative bond density wave (Relative BDW) . 
Although the operator pi^'_ (and therefore its expecta- 
tion value) does not vanish in the absense of the mag- 
netic field the relative BDW does not exist at / = (see 
Fig. 0) since the two conditions 34 < and (75 > can 
be resolved only when / > 0. 

All above long-range ordered states break sponta- 
neously translational symmetry of the underlying lattice 
(period doubling) and thus are doubly degenerate. Topo- 
logical excitations in these phases (^-kinks) carry unit 
charge Q — 1, as opposed to fractional charge Q — 1/2 in 
the quarter-filled case. This follows from the definition 
of the fcrmionic number carried by a single kink, 

/oo 
dx [Jr v {x) 4- J Lv (x)] 

u=a,p -°° 

dx d x (f>+(x) (15) 

T J-oc 

and the fact that each kink interpolates between the vac- 
uum values of the field <p + at x — > ±00 that differ by a 
period of the cosine potential, equal to wtt/2. 



D. Non-ordered phases 

In the previous subsection we have discussed the or- 
dered phases occuring under the conditions Eq. JTJJ, i.e. 
when both sectors in the effective Hamiltonian (|12fl ac- 
quire gaps. In all other cases there exist gapless excita- 
tions. These are characterized by correlation functions 
that at large distances decay as a power law 

(O(x)O(O)) ~ l/x d , 

where d is the scaling dimension of the operator O. Cor- 
relations with slowest decay are usually referred to as 
dominant. In the phase diagram Fig. [21 we categorize the 
gapless phases according to their dominant correlations, 
indicating the corresponding order parameter in paren- 
theses. In this subsection we briefly describe such phases. 

If the conditions Eq. I|14f) are reversed and K + > 1, 
K- < 1, then both sectors are gapless and the system is 
a Luttinger liquid. In this case the dominant correlation 
function is that of the pairing operator O sc at wavevector 
7r — 2kp\ 



(OUt,x)O sc (0)) 



cos[(7r — 2kp)x/a] 



There are two other cases when only one of the condi- 
tions Eq. I|14f) is violated. Then only one of the sectors 
acquires a gap while the other remains gapless: 

(i) If K + < 1, K- < 1, then the "charge" sector is 
gapped, but the "relative" sector remains gapless. By for- 
mal analogy with the Hubbard model we call this phase 
a Mott Insulator. In such state, incommensurate density 
or bond-density correlations with characteristic momen- 
tum 2kp 7^ 7r are dominant. Indeed, depending on the 
sign of c/4, either cos y/2~K(j) + or sin\/27r</>+ acquire finite 
expectation values. Therefore, either the transverse bond 
density p± or the 2k f part of the total charge density p+ 
display slowest algebraic decay of the corresponding cor- 
relation function determined by the "relative" sector (see 
Table [I] for bosonized expressions). So at (74 > 



(p ± (T,x)p ± (Q)) 



cos(2kFx/a) 



\K- 



(17) 



If 34 < 0, then Eq. (|17|) applies to the correlation function 
of p+. 

(ii) If K + > 1, JT_ > 1, then the "charge" sector is gap- 
less, but the relative sector acquires a gap. By analogy 
with spin-gap systems, we call such a phase a Luther- 
Emery liquid^ Now it is 6L that takes one of the two 
semiclassical values depending on the sign of 175 . It turns 
out, however, that this phase can only occur when g 5 < 0, 
so that {cosy/2~n9-) 7^ 0, and the dominant correlation 
is that of the pairing operator at zero momentum (with 
the power law determined by the "charge" sector) 



(Ol(r,x)O sc (0)) 



\V+T- 



1 1/*4 



(18) 



The phase boundaries as a function of U, V, /, r can be 
calculated by solving Eas. (|Al|l for when g 4 or g 5 is zero, 
or K + or K- is one. For completeness, these are writ- 
ten in Appendix lAl The complete weak-coupling phase 
diagram is plotted in Figure Eland was discussed in Sec- 
tion niBi 



IV+T — IX 



2/K 



+ \V-T — IX 



(16) 



IV. LOW ENERGY EFFECTIVE THEORY: 
CHAIN BASIS 

In this section we briefly review the effective low en- 
ergy theory that one can derive taking two independent 
chains as a starting point. Interchain hopping is then 
taken into account already at the bosonization level sim- 
ilarly to interaction terms. This approach is valid as long 
as t± <C vp(f)a, where vp(f) is the renormalized veloc- 
ity (see below). In the absence of the magnetic field the 
chain-basis description of the spinless ladder has been 
widely used in literature. 1516 Skipping inessential details, 
below we will give a brief review which will help to ana- 
lyze differences between the two approaches and further 
clarify the role of the magnetic flux. 

In the chain-basis approach, one starts by lineariz- 
ing the fermion dispersion on each chain in the vicin- 
ity of the two Fermi points, ±kp = ±7r/2, defines chiral 
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fermion fields, and then expresses the inter-chain hopping 
in terms of these fields. On the lattice, the magnetic field 
was introduced in the Hamiltonian (JJJ via the Peierls sub- 
stitution Eq. (J2J. Here it is convenient to split the phase 
exponential in Eq. (|2J) into its real and imaginary part. 
The real part contributes to the renormalization of the 
Fermi velocity, vp — > vp(f) = vp cos(7r/), which is of mi- 
nor importance as long as / is not too close to 1/2. The 
imaginary part can be written in terms of the densities 
of the left and right particles, so that the single-particle 
perturbation to free chiral fermions is of the form: 

Hi(x) = -h fl • J R -h L ■ J L (19) 

where the chiral densities (vector currents) are defined as 



Jfi = -R^a^Rp, 3 L = -tf a o a pLp. (20) 

Here a a (a = x, y, z) are the Pauli matrices and a and 
/3 are the chain indices. In Eq. I(19|) the vector currents 
appear to be coupled to the effective chiral "magnetic" 
fields 



h«= Qi±,0, -h\\), h L = (h±,0,hn 



with 



h\\ = At sin(7r/), h±_ = 2tj_. 
Bosonizing Eq. 1)19(1 directly one can find 



Hi 



■ cos v 2-k6_ sin v 2n9- 



too 



(22) 



Notice that only the relative field, (p- — (<f>i — (p2)/V^ 
and its dual, 6L = (8i — 9-i)j\f2 appear in Eq. l|2*2^l since 
we are discussing inter-chain processes. To avoid confu- 
sion, we remind that here (f>- and 9- are the differences 
between the corresponding bosonic fields defined at each 
chain. The remainder of the Hamiltonian in the relative 
sector originates from the kinetic term and the interac- 
tion, so the total Hamiltonian has the form 



H- = Kg + 7~Lsg + "Hi, 



(23) 



where Hg is the Gaussian model with the interaction 
parameter K expressed in terms of g s = {AV — U)a in 
the standard way, Eq. (|A4|) . and Hsg °c U cos is 
the sine-Gordon term. 

The last term in Eq. ((23(1 Hi is the sum of a nonlocal 
vertex operator and a derivative of the dual field . The 



latter appears only in the presence of the magnetic field 
and may be interpreted as the bosonized version of the 
Lorentz-invariant current j cx JJj — J£. It is important 
to realize, however, that this is not the physical relative 
current already discussed in Section IIII Al Moreover, j 
appears not to be gauge-invariant: its expectation value 
in the ground state of the non-interacting system is non- 
zero even in the absence of inter-chain hopping. In what 
follows we demonstrate that, physically, Hi [Eq. I|22|l ] 
describes the splitting of the two Fermi points in the 
chain basis into the four Fermi points introduced in Sec- 
tion nnsi 

Treating Eq. 1(23(1 presents certain difficulties related to 
the non-local nature of the perturbation 1(22(1 . Global ro- 
tations of the quantization axis for the Abelian bosoniza- 
tion, which proved to be efficient at / = (see Refs. Il6l 
IT^ l. arc not helpful here because of the more compli- 
cated structure of Hi . Let us therefore perform a chiral 
rotation which makes the "magnetic" fields h# and 
antiparallel and aligned along the z-axis: 



h' fl = (0,0,-fc), h / i = (0,0,ft), 



(24) 



where h 



hfi. This simplifies the bosonic form 



(21) Eq. 1(22(1 which now contains only the derivative term 



Hi 



-h\l-d x 9- 

IT 



(25) 



To achive this result we need to rotate the right field 
by u> and the left field by — lo about the y axis where 
tana; = h±/hu (so that ui = it/2 corresponds to the 
absence of flux). In terms of the fermion operators this 
chiral SU(2) rotation can be written as 



R a 



iu)o v /2 



a0 



Rp, L c 



— iuio y l1 



a0 



L . (26) 



The transformation 1(26(1 is closely related to the Bo- 
golyubov transformation Eq. (J3J) used to diagonalize the 
single-particle Hamiltonian in the band approach. Below 
we analyze the relation between the two in detail. 

While the rotation 1(26(1 simplifies the single-particle 
terms in the Hamiltonian, the interaction terms undergo 
a nontrivial modification. Since the rotation involves 
only "relative" fields, the "charge" sector is left unaf- 
fected and is a usual sine-Gordon model, with a cosine 
term cx — U cos v / 87r ( I > +. For the rest of this Section we 
will assume this interaction to be strong enough to gen- 
erate the gap in the "charge" sector and focus on the 
"relative" degrees of freedom. The Hamiltonian density 
describing the "relative" sector transforms to 
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n- = — 



K_U 2 _ + —{d x ^) 2 



■ cos v 87r0 



7T 3 / 2 Q!o L 



(7ra ) 2 (Trao) 2 
(0*0*) sin VSt:^ + (0 X <^ ) sin \/87n^ 



cos v 87r6> 



hJh x t 



(27) 



with all the coupling constants listed in the Appendix 1X1 
In Eq. (|27|l the relative sector represents a Z4 modelSLSSi 
modified by the g r term. Notice that the latter has zero 
conformal spin and scaling dimension 

d r = l + i(W^)>2 

for any K_ 7^ 1. Thus this term is irrelevant in the RG 
sense. Thus, as for the ladder in the absence of the mag- 
netic flux, the effective model appears to be two-cosine Z± 
model with a topological term. The only but important 
difference with the / — case is a nontrivial dependence 
of the coupling constants on the ratio t / shnrf through 
the rotation angle to. 

At h = 0, the Z4 model in Eq. Q27[) always displays a 
strong-coupling regime in the infrared limit accompanied 
by a dynamical generation of a mass gap. If K < 1, then 

term is relevant while the term gg with the dual field 
is irrelevant. Hence the field 0_ gets locked and the term 
hd x 0- has no effect on the corresponding long range or- 
der. In particular if g$ > 0, then the ground state is the 
relative CDW, which becomes the OAF when g^ changes 
sign (e.g. due to the variation of the flux) . However, lock- 
ing of the field 0_ does not prevent the gradient of the 
dual field d x 9_ to acquire a finite expectation value and, 
thus, produce a finite splitting of the Fermi momenta 
proportional to h (see the discussion below). 

In the opposite case K > 1, the gg term is relevant and 
leads to locking of the dual field 9_ . As long as h remains 
smaller that the gap generated in the "relative" sector, 
the ground state is the commensurate CDW, which re- 
places the MI shown in the phase diagram [3] In this 
phase the vacuum value of 9- resides in one of the min- 
ima of the cosine potential and remains spatially uniform, 
implying that (d x 9-) = 0. Thus no band splitting occurs 
in this regime reflecting the fact that the chains remain 
effectively decoupled. This is a manifestation of the An- 
derson confinement 1 — - the in-chain correlations suppress 
single-particle tunneling between the chains in the low- 
energy limit. Consequently the flux has little effect on 
this phase. 

The situation changes when h reaches a critical value 
h c proportional to the mass gap 

h c ~ A s /v F a ~ exp(-27rw F /|5r s |) < 1 (g s < 0). (28) 

At this point a commensurate- incommensurate transition 
takes place, the long-range order disapears and a finite 
gradient (d x 9 s ) emerges in the ground state, following the 



universal square-root increase ~ y/h— h c slightly above 
the threshold^ The commensurate CDW gets replaced 
by the MI phase with incommensurate leading correla- 
tions. The appearance of a finite average (d x 9 s ) indicates 
that the two doubly degenerate Fermi points, character- 
izing the bare-particle spectrum in the chain basis, are 
getting split. 

Note, that the commensurate-incommensurate transi- 
tion can be seen within the band basis approach if one 
takes into account an additional interaction process men- 
tioned previously in Eq. i|ll|) . which is usually disre- 
garded since formally it does not conserve momentum. 
The term is the interband backscattering 

R^RpL^pLa + h.c. (29) 

In bosonized form this can be written as 

, n 97 ,„ cos(V8^>„ + 2hx) (30) 
(27ra r 

where the relation between g-j and the microscopic pa- 
rameters of the model is given in Eq. |Alg| | . Adding 
this term to the bosonic two-band Hamiltonian i|12|) and 
making a shift 

(f)- -> (j>- - hx/V2^ (31) 

one transforms the relative part of the Hamiltonian 112|) 
to the form identical to the Z4 part in Eq. I|27|l (up to du- 
ality transformation 0_ <-» 6- ). The above analysis holds 
identically for this representation. While the case 
may be easier to describe within the chain approach, the 
band picture should always give the correct result. In 
particular, to find the above CDW order, one has to no- 
tice that for sufficiently small h when (cos ^/2tt(/)^ ) =^ 0, 
the 2k p component of the CDW order parameter listed 
in Table [I] acquires a nonzero expectation value (now 
2kp — 7r), and the MI phase in the bottom left of the 
phase diagram |31 becomes a CDW phase with long range 
order. 



V. STRONG COUPLING 

In this Section we discuss the behavior of the sys- 
tem in the strong-coupling limit, |J7|,|V| ^> t . In 
the atomic limit when hopping is completely neglected 
(<o = t± = 0), the particles are localized on sites, and 
there are four possible ground states. Which state has 
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the lowest energy is determined by the signs of the in- 
teraction parameters U and V: (i) if U, V > 0, then 
the ground state is the Relative CDW (the SDW Z in the 
"spin" language) depicted in Fig.^ (h) when U < and 
V > the state is the CDW also shown in Fig. ^ (hi) 
in the opposite case U > and V < all particles fully 
occupy one chain keeping the other empty; and finally, 
(iv) when U, V < we have complete phase separation. 

Of the above four ground states the first two are acces- 
sible in the weak coupling approach as can be seen from 
the phase diagram in Fig. [3J The phases (iii) and (iv) 
do not have the lowest energy when the bandwidth to is 
greater or of the same order as U and V and thus have 
no analog in weak coupling. 

To make further links with the weak coupling approach 
we now need to take into account the hopping terms and 
the magnetic field. Of the above four cases only the first 
one needs to be discussed in detail. Indeed, the last two 
do not appear in weak coupling, while in the case (ii) 
we have either doubly occupied or empty rungs, so that 
inter-chain hopping and the flux do not afect the proper- 
ties of the ground state. This ground state (CDW) was 
discussed in the previous section (this is the case where 
h is smaller than the gap, see Eq. (|28|l . 

Consider the limit where U > is the largest scale in 
the problem. Then we can project out states with doubly 
occupied rungs. Then at half-filling and at energies well 
below the local charge gap, there only remain configura- 
tions with exactly one fermion per rung. Accordingly, the 
relative degrees of freedom can be conveniently described 
in terms of local spin- 1/2 variables using the correspon- 
dence between two single-fermion states at a given rung 
n and the eigenstates | f) and | J.) of the operator S*. 
The standard Schrieffer- Wolff transformation 2 — leads to 
the following effective spin-chain model: 



Keff = JoJ2[l( e2iwfS nS- + , + e- 2i ^S-S+ +1 ) 



+AS^S^ +1 + h±S* 



(32) 



where Jo = 2£q/L/ is the exchange constant, A = (Jo + 
V)/Jq is the anisotropy parameter, h± = t±/Jo is the 
transverse field and / is the external flux in the origi- 
nal ladder model. Notice, that the Schrieffer- Wolff tran- 
sormation leaves the "charge" sector gapped, so that all 
subsequent analysis pertains to the "relative" sector. 

In the spin language, the Relative CDW order param- 
eter (see Table QJ corresponds to the staggered magneti- 
zation in the z-direction (— l) n S*, while the OAF order 
parameter is the staggered magnetization in the y direc- 
tion, (— l) n SV. Magnetization in the x direction corre- 
sponds to the transverse bond density. Since the uniform 
transverse field hx breaks the U(l) symmetry of the XXZ 
chain, the staggered component (— l) n S£j never acquires 
a non-zero expectation value. 

The effective spin model Eq. l|32[l is not integrable 
and its general solution remains unknown. Nevertheless, 
there exist at least three cases where further progress 



can be made: (a) the case / = which has been stud- 
ied previously (see e.g. Ref . |30S) : (b) the vicinity of the 
S't/(2)-symmetric point A = 1 and / < 1; and (c) the 
vicinity of / = 1/2. In what follows we discuss these 
three cases. 

(a) In the absence of the flux / = the Hamilto- 
nian (|32|l becomes equivalent to the transverse field XXZ 
model. Despite not being integrable, much is known 
about such models^ Here, we summarize the results for 
the sake of completeness referring the interested reader 
to the literature for more details. 

First consider the case V > (so that A > 1) and 
the limit V ^> 2t^/U in which the exchange anisotropy 
becomes very large, A ^> 1. In this case one can retain 
only the last two terms in (|32() . so that the Hamiltonian 
becomes equivalent to the one-dimensional Ising model 
in a transverse magnetic field. At h± < A/2 the ground 
state is the ordered Neel phase with ((-l) n S*) ^ 0. This 
ordering translates to the Relative CDW for the origi- 
nal ladder, in agreement with the weak-coupling picture. 
However, when h± = A/2 a Z2 (Ising) transition to a dis- 
ordered phase takes place. This transition is not present 
in the weak-coupling phase diagram, so we will not dis- 
cuss it any further. 

Now, if V < but is sufficiently small (-1 < A < 1), 
then in the absence of h± the Hamiltonian (|32|l corre- 
sponds to the critical XXZ model, which has a gapless 
excitation spectrum and displays dominant antiferromag- 
netic fluctuations in the xy-plane. A non-zero transverse 
field breaks the U(l) symmetry down to Z2 partially po- 
larizing the spins in the x-direction, so that ((— l) n S%) 
develops an expectation value. This is the OAF phase 
already discussed in connection with the weak-coupling 
phase diagram. 

Finally, if A < —1 then the model is the easy axis XXZ 
ferromagnet with (S*) ^ 0. Such a ground state corre- 
sponds to the case (iii) above, i.e. all particles localized 
on one of the chains. 

(b) We now consider the case |A — 1| <C 1 and /< 1. 
In this limit the Hamiltonian i|32[l can be represented (to 
lowest order in /) in the form 

H = Jo 22 {^™ ' &n+i + (A — tySn^n+i 

n 

+ h {l (S n x S n+1 ) z + h ± S*}, (33) 

i.e. we have a model of a weakly anisotropic XXZ spin 
chain in a magnetic field h± = t±/Jo along the x direc- 
tion, and also perturbed by the term proportional to the 
z-component of the spin-current (S n x S n+ i) z with an 
amplitude proportional to the flux, h\\ — irf. Bosoniz- 
ing around the S'f7(2)-symmetric point in the standard 
wayji& one can show that the last two terms in Eq. I)33|) 
can be written in the form Eq. (|19fl with the chiral "mag- 
netic" fields, Eq. 121|) . rcnormalized by the new "band- 
width" Jo. Thus we find that the bosonized form of 
Eq. (|33|) has the structure of Eq. (|23|l . where the first 
two terms Tic an< i 7~Isg constitute the Abelian bosonic 
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representation of the SU (2)-symmetric Heisenberg model 
in the scaling limit, with the Lutting parameter and 
the effective coupling constant renormalized by the weak 
anisotropy term |A — 1| cx V in Eq. I|33l) . 

Notice that the resulting Hamiltonian is basically the 
same as that obtained at weak coupling in the chain ba- 
sis (|23|) . Therefore the next step, namely the rotation 
Eq. d2HJ and its result Eq. E3 can be performed in ex- 
actly the same manner as before. As a result, at least in 
the region U > 0, all conclusions drawn at weak coupling 
also hold true in the strong coupling limit. In particular 
the region of the phase diagram discussed in the previ- 
ous Section survivies (up to renormalizations of the phase 
boundaries) in the strong coupling regime as well. 

(c) Finally, we consider the case where flux is close 
to one half. Here we are close to the transition to a 
band insulator. This region of the phase diagram can 
not be well treated in our weak coupling approach. In 
this sense, the strong coupling arguments compliment the 
weak coupling picture presented in Section IlIII 

Consider a gauge transformation 



C+ __. C+„2i7r/n q- c+„-2i7r/n 

which transforms the Hamiltonian (|32|) to 



(34) 



+h ± (S+e 2l7rfn + S-e- 2mfn ) 



(35) 



The model (|3*5|) is completely equivalent to Eq. but 
now the transverse field is non-uniform. The situation 
significantly simplifies when \ f — 1/2| <C 1. In the case 
this field becomes almost staggered and can be directly 
bosonized (in the region |A| < 1, i.e V < 0).- 6 As a 
result, Eq. I|35|l becomes 



H 



dx 



ATI 



1 

K 



(36) 



dx cos 



y/nQ(x) - 2vr 



/ 



where the Luttinger liquid parameter is given by 

K = n/2(n — cos -1 A). (37) 

For / = 1/2, the cosine term is relevant and generates 
a gap in the spectrum. The resulting ground state is the 
band insulator and the gap corresponds to that seen al- 
ready in the single-particle problem in Section l!!! Al This 
can be seen from the fact that cos(y / 7r0) (which gains a 
non-zero expectation value at / = 1/2) corresponds to 
the uniform bond density. 

As / decreases away from 1/2, eventually the gap 
will close via a commensurate-incommensurate transi- 
tion. The system will now have gapless excitations in 
the relative sector (i.e as in a MI), in qualitative agree- 
ment with the weak coupling phase diagram Fig|3J 



Finally, if V > so that A > 1 one can still bosonize 
the Hamiltonian (|35|l . although now there is an extra 
term proportional to cos \/167r$ which is relevant (for- 
mally, this is the case K < 1/2 which is not captured 
Eq. (|37|) ). We now have two competing cosine terms in 
the Hamiltonian, cos Vl67r$ and cos[\/nQ(x) — 27r(l/2 — 
f)x/a]. The dual field term has the smallest scaling di- 
mension and therefore, when / = 1/2, determines the 
character of the ground state. If the flux is decreased, 
then again the order is destroyed. In this case, however, 
the other relevant operator in the problem, cos VT&tt$ 
should at this point acquire a non-zero expectation value. 
This is the transition between the band-insulator and the 
relative CDW, again in agreement with Fig El 



VI. SUMMARY 

We have investigated a model of interacting spinless 
fcrmions hopping on a two-leg ladder in the presence of 
an external magnetic field at half-filling. Using bosoniza- 
tion techniques, we constructed the effective low energy 
theory where the coupling constants acquired non-trivial 
dependence on the external flux. Consequently the flux 
results in several phase transitions shown in the weak- 
coupling phase diagram Fig. [3] i.e. BKT transitions be- 
tween ordered and disordered phases and U(l) transi- 
tions between different ordered phases. 

Furthermore, we extended our weak-coupling picture 
by the special consideration of the case of weakly cou- 
pled chains at small flux. We solved the corresponding 
effective theory using the chiral rotation Eq. I|26|) . As a 
result we described the commensurate-incommensurate 
transition from the MI phase shown in Fig. [21 at U < to 
the CDW. This transition happens when the parameter 
\J t 1 + sin 2 7r/ becomes small enough [see Eq. d28fl ]. 

The weak coupling analysis is complimented by the 
strong coupling arguments. In particular we showed 
that in some cases (in particular, for small flux and 
small in-chain interaction V) , that the weak-coupling ap- 
proach and strong-coupling approach lead to identical 
low-energy theories. 

Finally, we discussed the persistent (diamagnetic) cur- 
rent flowing in the ladder when the external magnetic 
field is applied. We showed that even in the limit of 
small / and small t±, all electrons in the system con- 
tribute to this current and therefore it is a non-universal 
feature which can not be described in the traditional 
field-theoretic approach. We speculate that it may be 
an interesting physical quantity to investigate within the 
non-linear bosonization scheme. 31 
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APPENDIX A: RELATION BETWEEN 
PARAMETERS AND PHASE BOUNDARIES 

Here we list the effective coupling constants for the 
effective low energy theory and their relation to the mi- 
croscopic parameters in the original ladder Hamiltonian 
Eq. HJ. 

First, we show the parameters of the g— ology for the 
effective interaction Eqs. (|l(Jfl and 111|) in Section lill Bl 



.91 



COS 2 71"/ — T 2 

sin 2 7r/ + t 2 



Ua tan 2 nf + 2Va 



COS 2 71"/ 



g 2 = Ua + 2Va 



COS 2 71"/ 



.93 = 2Va + - 



1 



sin 7r/ + t 2 



Uar 2 

COS 2 71"/ 



(Ala) 



(Alb) 



(Ale) 



+2Va(sin 2 nf - t 2 tan 2 tt/) 



.94 = 52; 



.95 = .9i; 



4 \ sin 7rj + 



.97 



COS 2 71"/ 



2Va 



sin 2 7r/ + r 2 



(Aid) 



(Ale) 



(Alf) 



(Alg) 



The two "Luttinger parameters" are then given by 



K+ = 



l-(gi± g 3 )/47 



1 + (51 ± <? 3 )/47 



(A2) 



The phase boundaries in Fig. (within the accuracy 
of one-loop renormalization group approach) are given 
by the following conditions: 



,9i + .93 = 
,9i - .93 = 
.94-0 
35 = 



U 
2V 



2 cos 2 nf - t 2 



COS 2 71"/ 

U _ sin 2 tt/(2 cos 2 tt/ - t 2 ) + r 4 
2T^ sin 2 7t/ (cos 2 7t/ — r 2 ) — r 2 

2V 
U 
2V 



cos 2 7r/ 

r 2 

• 2 <• 
sin 7T/ 



(A3) 



Notice that the boundaries depend only on the ratio 
U/V. 

Finally, we list the effective constants for the effective 
low energy theory in the "chain basis" Eq. I|27|) 



[2Va sin 2 oj + Ua cos 2 u] ; 



2Va - Ua 



sin 2 ui: 



g r = (2Va — Ua) sin2w; 
g s = [4Va cos 2 u> + Ua (sin 2 uj — cos 2 w)l , 
1/2 



1 - gs/2-KVs 



1 + g s /2irv s 



(A4) 



APPENDIX B: SINGLE PARTICLE PROPERTIES 

Here we point out trivial properties of the single- 
particle spectrum © for the sake of completeness. The 
upper band has its absolute minimum at k — with 

{£p)min = ep{0) = T - COS 71"/. (Bl) 

The lower band has its absolute maximum at k — ±n 

(e a )max = e Q (7r) = cos7r/ - t. (B2) 

Thus, if cos 7r/ < r, the spectrum exhibits a single- 
particle gap and the non-interacting system is a band 
insulator at half filling. 
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FIG. 6: The phase diagram of the non-interacting system. 
The vertical axis is the external flux and the horizontal axis 
is the ration between inter-chain and in-chain hopping am- 
plitudes. Fig. 0]demonnstrates the four typical shapes of the 
spectrum corresponding to the four corners in which the phase 
diagram is separated by the two lines Eq. 1B4II . 

If the flux is not too small then there is a no n-trivial 
extremal point cosfci = — cot ixj ' \/ t 2 + sin 2 7r/. For the 
upper band the energy e^(fci) has to be compared to the 
remaining extremum 6/3(71"). If sin 2 7r/ > tcos7t/ the 
band acquires a double-well shape as shown in plots II 
and IV in Fig.0|and k\ turns out to be the absolute max- 
imum for the uper branch of the spectrum, while k = ±7r 
are local minima. For the lower branch of the spectrum 
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and under the same condition, k\ is the absolute mini- 
mum, while k — is the local maximum with the energy 
value 

e«(0) = -t-costt/ (B3) 

The two aforementioned conditions, namely, 

cos7r/2 = t; sin 2 7r/i = tcostt/i (B4) 

define the two boundaries in the phase diagram of the 
non-interacting system shown in Fig. These two lines 
separate the phase diagram of the non-interacting system 
(Fig. into four parts, where the spectrum has one of 
the four shapes shown in Fig. 0] 

APPENDIX C: BOSONIZATION CONVENTIONS 

Here for completeness we define our bosonization con- 
ventions. The chiral bosonic fields are introduced via the 
correspondence 

(R„{x),L v {x)) -> (k„/V2^) e*^+? h \ (CI) 
Jrv = d x cj> Rl ,/sJ~n, J Lu = d x ^L V j\pK. (C2) 



Here v — a, (3, ao is an ultraviolet cutoff in the bosonic 
theory, k v are Klein factors that ensure proper anticom- 
mutation relations between the fermionic fields with dif- 
ferent band indices in representation (|C1(I . The n u satisfy 

{ k pj k i>} — 25^^; = —%. (C3) 

In addition, we impose a nontrivial commutation relation 
between the right and left bosonic fields belonging to the 
same band: 

=<W*- (C4) 

The left- and right-moving fields can be combined into 
the field <j> and its dual counterpart 9 

<Pa = 4>a + <^«; ^a = 0a-^a: ( C5 ) 

with II = —d x 6 u being the momentum conjugate to <\> v . 
The linear combinations 

cb ± = (0 Q ± M/y/2; 9± = (6 a ±6 )/V2. 

describe collective bosonic degrees of freedom and the 
symmetric and antisymmetric sectors of the theory. 
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